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ABSTRACT 

The early stages of decelerating gamma-ray burst afterglow jets have been notoriously difficult to resolve nu- 
merically using two dimensional hydrodynamical simulations even at very high-resolution, due to the extreme 
thinness of the blast wave and high outflow Lorentz factors. However, these resolution issues can be avoided 
by performing the simulations in a boosted frame, which makes it possible to calculate afterglow light curves 
from numerically computed flows in sufficient detail to accurately quantify the shape of the jet break and the 
post-break steepening of the light curve. Here, we study afterglow jet breaks for jets with opening angles of 
0.05, 0.1 and 0.2 radians decelerating in a surrounding medium of constant density, observed at various angles 
ranging from on-axis to the edge of the jet. A single set of scale-invariant functions describing the time evo- 
lution of afterglow synchrotron spectral break frequencies and peak flux, depending only on jet opening angle 
and observer angle, are all that is needed to reconstruct light curves for arbitrary explosion energy, circumburst 
density and synchrotron particle distribution power law slope p. These functions are presented in the paper. 
Their time evolutions change directly following the jet break, although an earlier reported temporary post-break 
steepening of the cooling break is found to have been resolution-induced. We compare synthetic light curves 
to fit functions using sharp power law breaks as well as smooth power law transitions. We confirm our earlier 
finding that the measured jet break time is very sensitive to the angle of the observer and can be postponed 
significantly. We find that the difference in temporal indices across the jet break is larger than theoretically an- 
ticipated and is about -(0.5 + 0.5/?) below the cooling break and about -(0.25 + 0.5/?) above the cooling break, 
both leading to post-break slopes of roughly about 0.25 — l.3p, although different observer angles, jet opening 
angles and heuristic descriptions of the break introduce a wide range of temporal indices. Nevertheless, the 
post-break slope from our constant density ISM simulations is sufficiently steep to be hard to reconcile with 
post-break slopes measured for the Swift sample, suggesting that Swift GRBs mostly do not explode in a ho- 
mogeneous medium or that the jet breaks are hidden from view by additional physics such as prolonged energy 
injection or viewing angle effects. A comparison between different smooth power law fit functions shows that 
although smooth power law transitions of the type introduced by Harrison et al. 1999 often provide better fits, 
smooth power law transitions of the type introduced by Beuermann et al. 1999 or even sharp power law fits are 
easier to interpret in terms of the underlying model. Light curves and spectral break and peak flux evolution 
functions will be made publicly available on-line at |http : / /cosmo . nyu . edu/af terglowlibrary| 

Subject headings: gamma-ray bursts: general - hydrodynamics - methods: numerical - methods: data analysis 
- shock waves 



1. INTRODUCTION 

Gamma-ray burst (GRB) Afterglows are produced by non- 
thermal radiation from collimated decelerating relativistic 
outflows following the collapse of a massive star or a neu- 
tron star-neutron star or neutron star-black hole merger (for 
reviews, see e.g. |plraiill2004t iMeszarosll^M lN^al2007l 
lGranotll2"007l) . Because they originate from cosmological dis- 
tances, their jet nature can not be observed directly but is ex- 
pected theoretically from constraints on the energy budget of 
the outflow (the isotropic equivalent energy of the afterglow 
often being comparable to the solar rest mass) and inferred 
observationally from the jet break in the light curve. This 
break has been observed at various wavelengths ranging from 
radio to X-rays, and marks the onset of a steepening of the 
decay of the light curve. In this paper we present the most 
accurate description to date of the temporal and spectral evo- 
lution of the afterglow signal during the jet break, based on 
detailed relativistic hydrodynamics (RHD) calculations of the 
afterglow blast wave decelerating in a homogeneous medium. 

The steeper decay following the break is the result of two 



changes in the outflow. On the one hand the jet is start- 
ing to become less collimated. As a result, the area of the 
blast wave front increases and the jet decelerates faster than 
before because it starts to sweep up more circumburst mat- 
ter. On the other hand, the ongoing deceleration even with- 
out spreading reaches a point where the relativistic beaming 
cone of the synchrotron emission at and behind the shock 
front becomes sufficiently wide for the lack of flux beyond 
the edges of the jet to become visible, whereas before only 
a small patch along the direction of the observer could be 
observed and a jet was indistinguishable from spherical out- 
flow. Both effects are expected to occur approximately around 
the same point in time, when jet half-opening angle 6q ~ 
I/7, with 7 the fluid Lorentz factor behind the shock. For 
the sideways spreading this is because #0^1/7 marks the 
point where the fast spreading in the frame comoving with 
the jet becomes noticeable as well in the frame of the ob- 
server, while for the edge effect it marks the point where the 
beaming cones have become as wide as the jet itself. Al- 
though the wideni ng of the jet was originally argued to be 
the stronger effect (Rhoa dsTl999h . subsequent numerical stud- 
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ies dGranot et al.ll200UlKumar & Granolll2003t iMeliani et all 
20071: IZhang & MacFadvenl l2009t lVan Eerten et alj 1201 lbt 
Wvgoda et alj|201 It iDe Colle et al.ll2012al) reveal the spread- 
ing of the jet to be more steady than the runaway process 
from iRhoadsl (119991) . for observationally relevant jet open- 
ing angles. As a result, the edge effect plays a strong role 
in shaping the jet break and the angle of the observer rela- 
tive to the jet axis becomes relevant even for observer angles 
within 9n (I Van Eerten et aTll2010at IVan Eerten & M acFadven 
l2012ttlVanEertenetal]l2011ah . Bv now, theoretical mod- 
els incorporate more realistic descriptions of jet spreading 
dGranot &Piranll2012l) . 

Observational signatures of collimated afterglow outflow 
for a number of GRBs first started to emerge in the late 
nineties, both from the overall steepness of the light curve 
compa red to theoretic al expectations for a spherical explo- 
sion dSari et alJ 119991) and observations of the jet break 
dBeuermann et al.l 1 1 999fc IFruchter et alJl!999t [Harrison et al. 
1999; Kulkarni et alj 1 1 999b . Because of the complexity 
of the dynamics of decollimating jets, afterglow jet breaks 
have been modeled by heuristic functions for the purpose 
of data fitting. From the beginning c onnected power laws 
have been used by many groups (e.g IFruchter et al.1 119991 : 
iKulkarni et alJI 1999f) but also power laws with a sm ooth tran- 
sition (Beuermann et al. 1999; Harrison et alj!999l using dif- 
ferent descriptions for the transtion). Sharp and smooth 
power laws to describe jet breaks (or breaks in general) 
in afterglows have also been used in many more recent 
studies (e.g. IZeh et al J 120061: iLiang et all 120091: lEvans et all 
2009t IRacusin et al] 12009b iNicuesa Guelbenzu et alJ 1201 It 
Oates et al.l 1201 It iFong et al] 120121; iPanaitescu & Vestrandl 
20121) . The advantage of a general heuristic function to de- 
scribe afterglow breaks is that they do not necessarily as- 
sume an underlying model (i.e. jet break) but aim to de- 
scribe the observed shape of the data in a concise and con- 
venient manner. Identifying breaks as jet breaks is a sep- 
arate step, where the steepening of the break and pre- and 
post-break relations between s pectral and temporal slope (the 
"closure relations", see e.g. Zhang & Meszlroi 120041) are 
compared against theoretical expectations. Using such meth- 
ods a lack o f afterglow breaks has been reported for the 
Swift sample dKocevski & Butlerl 120081: IRacusin et alj 120081 
2009), which has b een attributed to the quality of the data 
(Curr anet alj |2008[) and the effect of the observer angle 
( Van Eerten et alj 1201 lal providing a mechanism by which 
jet breaks can be postponed b eyond Swift's capability to ob- 
serve, as suggested by e.g. iKocevski & Butlerl 120081) . Re- 
cent light curves from numerica l simulations demonstrate 
dVan Eerten & M acFadven 2012a) that the shape of the after- 
glow synchrotron spectrum changes strongly directly follow- 
ing the jet break, which renders the standard application of the 
closure relations unreliable and might serve to e xplain the lack 
of suc ces in using them to identify jet breaks dRacusin et alJ 
I2009T) . 

It ha s recently been pointed out (IVan Eerten & MacFadvenl 
12012a!) that the shape of the jet break in the light curve is de- 
termined by a scale-invariant function that depends only on 
initial jet opening angle 8q and observer angle 9 t, s and that 
scales in a straightforward manner between jet energies and 
between circumburst densities. This function is calculated 
from high-resolution relativistic hydrodynamics (RHD) simu- 
lations of the jet dynamics in 2D that include lateral spreading 
and deceleration to trans-relativistic velocities. 

This scale invariance has a number of useful practical im- 



plications. It makes it possible to distill a description of the jet 
break from numerical simulations that includes the full com- 
plexities of 2D trans-relativistic jet dynamics without the need 
to explicitly probe the parameter space in burst explosion en- 
ergy and circumburst density with time-consuming RHD sim- 
ulations. The resulting jet break description will be general 
and uniquely constrains the post-break closure relations and 
light curve slope. Existing smooth power law descriptions 
of the break can be compared against the simulation-derived 
shape. Simulation-derived jet break functions can even be fit- 
ted directly against the data in order to identify jet breaks. 
When these dimensionless jet break functions are scaled in 
order to fit the real time evolution of the data, the ratio po /Ej SO 
is obtained, yielding important constraints on the physics of 
the progenitor. 

Although significant progress has been made recently 
in numerically resolving afterglow jets properly using 
adapt i ve mesh refinement techniques dZhang & M acFadven 
[20091: IVan Eerten et al l EOlOal 1201 lbt I Wveoda et al.1 1201 It 
IVan Eerten et al.l 120 12t IDe Colle et al.l l2012al) . to date the 
early stages of the blast wave evolution have not been 
fully resolved due to the extreme sharpness of the blast 
wave profile in the sel f-similar Blandford-McKee (BM, 
lBlandford~& McKee 119761) solution for an ultra-relativistic 
blast wave that provides the initial conditions for the sim- 
ulations. As a consequence of this steepness, most matter 
in the blast wave is contained within a thin shell of typical 
width AR ~ R/I2j 2 , with R the blast wave radius. When this 
thin shell is not completely resolved, this leads to a transient 
startup phase characterized by a temporary artificial drop in 
Lorentz factor of the outflow. Only after the blast wave has 
been evolved for some time does the fluid profile return to the 
shape predicted from analytically evolving the initial condi- 
tions. In practice, this transient feature would typically still 
be present at least at the onset of the jet break, because of the 
trade-off between decreased resolution at earlier starting times 
(due to the 7-dependency of the width) and decreased validity 
at late times (the closer to the jet break, the less valid the as- 
sumption of purely radial flow). In addition, when computing 
synthetic light curves, one has to account for the fact that the 
observed flux at a given point in observer time is made up of 
emission from a wide range of emission times, with contribu- 
tions from the back of the jet being emitted earlier than those 
from the front in order to arrive at the same time. 

In order to obtain a truly accurate picture of the shape of the 
jet break, it is therefore required to resolve the pre-break flow 
of the blast wave completely up to sufficiently high Lorentz 
factor that the effect of the start-up transient is removed. In 
the current study we completely resolve the afterglow blast 
wave at extremely high Lorentz factors and early times by 
performing the RHD calculation in a different frame than the 
usual burster frame, which is at rest with respect to the ex- 
plosion engine and the observer (aside from a cosmological 
redshift). By changing to a frame moving at fixed relativistic 
velocity along the jet axis, the narrowness of the jet profile 
due to Lorentz contrac tion is reduced and all relative L orentz 
factors become small (MacFa dven & van Eerten|[2012l) . The 
price that is paid for this frame transformation, the loss of 
simultaneity across the grid, can be accounted for when the 
radiation from the evolving blast wave is calculated. 

The features of the dynamics of narrow jets and ultra-high 
initial Lore ntz factor (> 100) flows will be p resented in a sep- 
arate study (MacFad ven & van Eert en 2012). In this work we 
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limit ourselves to the radiation from afterglow jets and de- 
termine the general shape of the jet break for afterglow blast 
waves that start out highly relativistic (Lorentz factor of 100) 
and have an initial half-opening angle 6q of 0.05, 0.1 or 0.2 
rad., moving into a homogeneous environment. The observer 
angle is varied from observers looking straight into the jet to 
observers positioned on the edge of the jet. In §|2]we discuss 
the methods of our RHD simulations and our implementa- 
tion of the BM initial conditions in a boosted frame. In $3] 
we discuss how light curves are calculated from simulations. 
In ^4] we show our results for the small set of key character- 
istic quantities (i.e. the break frequencies of the power law 
synchrotron spectrum and the peak flux) that determine the 
afterglow spectrum. We then use the characteristic quanti- 
ties to calculate afterglow light curves at optical and X-ray 
frequencies in §0 and compare the shape of the jet break to 
earlier parametrizations from the literature. Our results are 
summarized and discussed in §|6] Some technical aspects con- 
cerning radiative transfer from a Lorentz -boosted simulation 
frame are discussed in appendixlAl 

2. METHODS FOR BLAST WAVE DYNAMICS 

We assume that the radiation and the dynamics of the colli- 
mated relativistic blast wave can be separated. This assump- 
tion remains valid as long as the emitted energy is only a 
small fraction of the blast wave energy and as long as there 
is neglible feedback from the radiation on the jet dynamics. 
Additionally, we assume that the magnetic fields generated at 
the front of the blast wave also contain only a small fraction 
of the available energy. Under these assumptions the jet dy- 
namics can be computed using RHD simulations. 

2. 1 . Description of the RHD code 

We employ the RAM parallel a daptive-mesh refinement 
(AMR) code (IZhang & MacFadyenll2006h . The AMR tech- 
nique, where the resolution of the grid can be dynamically 
doubled locally where necessary, is important in order to re- 
solve the wide range of spatial scales involved, given the 
kR~R/\2Y 2 width of the blast wave in the lab frame in 
which the origin of the explosion is at rest, where 7 can be 
> 100 for a typical afterglow b last wave. Ram make s use 
of the PARAM ESH AMR tools (fMacNeice et afl l2000h from 
FLASH 2.3 dFryxell et al.1 120001). We use the second order 
F-PLM scheme (Zhang & MacFadven 2006) for the hydrody- 
namical evolution. In this study, a Taub equation of state is 
used where the adiabatic index smoothly varies between 4/3 
for relativistic fluids and 5/3 for non-relativistic fluids, as a 
function of the ratio between comoving density a nd pressure 
(Mi gnone et alJ[2005HZhang & MacFadvenll2009h : 

(h-4p)(h-p)=p\ (1) 

where p the pressure, p the comoving density and enthalpy 
h = pc 2 + p + e, with e the energy density. 

2.2. Scale invariant initial conditions 

Blast waves for three different jet half opening angles were 
simulated for this study: #0 = 0.05, 0.1 and 0.2 rad. The cir- 
cumburst number density no of the interstellar medium (ISM) 
is kept fixed at 1 cm" 3 , or po = m p g cm -3 , with density p and 
number density n related according to p = n x m p , where m p 
the proton mass. A more general expression for the circum- 
burst density environment is given byp ( > = po, re f(r/r re f)~ k 



A r , with r the radial coordinate, po,re/> r re f and A parame- 
ters setting the density scale and k setting the power law slope 
of the medium. Boosted frame simulations of blast waves de- 
celerating in a stellar wind environment where k = 2 will be 
presented in a follow-up study. All simulations start from the 
self-similar BM solution for impulsive energy injection with 
isotropic equivalent explosion energy Ei so set at 10 53 erg. The 
actual values for the initial energy and circumburst density are 
completely arbitrary and the hydrodynamics equations can be 
expressed in terms of dimensionle ss variables. Generaliz ing 
these variables from the ISM case (Van Eerten et al. 2 0121) to 
arbitrary k values, we have 



that are scale invariant under the transformations 

F' — kF- 

A' = XA, 

r' = (K/X) x ^- k \, 

t' = ( K /\f^- k) t. (3) 

All scale-invariance relations follow from straightforward di- 
mensional analysis, and are therefore not limited to the ultra- 
relativistic self-similar BM solution but apply throughout the 
evolution of the blast wave when jet spreading and decelera- 
tion occur. 

2.3. Simulations in a boosted frame 

Two challenging aspects of numerically simulating BM 
type outflows are the severe steepness of the radial profile 
of the various fluid quantities (i.e. the primitive quantities 
Lorentz factor 7, comoving density p, pressure p, and con- 
sequently the conserved quantities as well) and the ultra- 
relativistic nature of the outflow. Resolution issues regarding 
numerically calculated blast wav es and light curves have been 
discussed by various authors (IZhang & MacFa dven! 120091 : 
iVan Eerten etaTl 120 10U I20TI [Pe Colle et al.1 12012al) . As 
mentioned in the introduction, the most striking feature of an 
underresolved BM blast wave is a temporary spurious drop 
in Lorentz factor near the shock front. Because the observed 
flux strongly depends on the Lorentz factor (F„ oc ■y 2 , due to 
relativistic beaming), this strongly impacts the light curve. In 
order to understand the early time dynamics, it is important to 
start from a time when outflow peak Lorentz factor 7 ^> 1/6q 
(the point at which sideways spreading is expected to become 
relevant and when the edges of the jet become observable) 
and ideally, any transient behavior due to numerical resolu- 
tion should have subsided before this point. 

In the current work we have used cylindrical (R, z) coordi- 
nates . The initial conditions we re provided by the BM solu- 
tion dBlandford & McKee 1976), but expressed in a Lorentz 
boosted frame (MacFadv en & van E erten 2012). For all sim- 
ulations in this study, the simulation frame was boosted with 
Lorentz factor 75 = 5. All jets start with peak lab frame 
Lorentz factor 70 of 100 at the on-axis tip of the jet, though 
some were also run with 70 = 50 to check for convergence 
(which is expected for 70 3> 1/ Oo)- 

2.4. Resolution 

The simulation frame time duration of each simulation was 
2 x 10 7 seconds. The cylindrical grids run from -2 x 10 7 Is 
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FIG. 1. — Comparison of jet (8q = 0.05 rad) resolution between BOXFIT 
(dashed curve) and moving frame (solid curve) for a jet moving into the ISM. 
Top figure shows the evolution of the peak blast wave Lorentz factor (using 
f} X 7, the spatial component of the four-velocity) along the jet axis. Accord- 
ing to the BM solution, 7 oc f~ 3 / 2 and this slope is indicated by a dotted line. 
The bottom plot shows a light curve with jet break observed in between v m 
and v c . 

(lightseconds) to 2 x 10 7 Is in the z direction and out to 2 x 10 7 
Is in the R direction perpendicular to the jet axis. The initial 
peak refinement level is 15, with 8 base level blocks and 8 
cells per block in each direction. The smallest cell size at peak 
refinement level is therefore 5z = 25R = 19.1 Is = 5.72 x 10 11 
cm. Note that these are expressed in the boosted frame, so 
that the resolution in the z-direction is better compared to the 
lab frameQby a factor of 75. We enforce an upper limit on the 
total number of blocks on the grid. As the blast wave expands 
in size on the grid, the peak refinement level is decreased in 
order not to exceed this block limit. 

The top panel of Fig. Q] shows the evolution of the Lorentz 
factor at the shock front (in the lab frame) along the jet axis, 
compared to that in earlier work (Van Eerten et al. 20 12|), for 
the case where Qq = 0.05 rad. The Lorentz factor is measured 
at the numerically determined momentum maximum, which 
is slightly behind the exact position of the shock front. The 
Lorentz factor of the boosted frame simulation (solid line) 

1 Throughout this paper we will use "lab frame" to refer to the frame in 
which the origin of the explosion and the unperturbed interstellar medium are 
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FIG. 2. — Early time (observer time t = 0.2 days) pre-break spectrum. Et sa = 
10 53 erg, n = 1 cnT 3 (ISM), p = 2.5, 6,. = 0.1, e B = 10~ 3 , £v = 1- Dashed grey 
lines indicate asymptotic slopes, from left to right: 2, 1 /3, (1 —p)/2, —p/2. 

agrees with the BM solution at that position to within ~ 1%. 
The dashed line shows the BM scaling 7 oc f~ 3,/2 appropriate 
for the behavior of the shock Lorentz factor. 

3. RADIATION 

The algorithm used to calculate the radiation for a given 
observer time, fre quency, angle and distan ce is nearly identi- 
cal to th at used in[\ ^n Eerten et al. (201 2|), which in turn was 
based on lSari et alj dl998l) and lGranot et alj dl999l) . The only 
difference is that it is now applied to a boosted frame simu- 
lation rather than a non-moving frame simulation. The radia- 
tive transfer equations are solved for a large number of rays 
through the evolving blast wave. The stepsize along the rays 
is set by the number of data dumps from the simulation (3000, 
although we found in practice that the light curves were con- 
verged even for 300 data dumps). The conceptual details of 
the radiative transfer approach for a moving frame are pro- 
vided in appendixlAl 

The emission and absorption coefficients are calculated for 
synchrotron emission and synchrotron self-absorption (s.s.a.). 
The local synchrotron emission spectrum is given by a series 
of sharply connected power laws, with peak flux and spectral 
breaks determined by the local state of the fluid. A relativistic 
distribution of shock-accelerated particles is assumed, carry- 
ing a fraction e e of the local internal energy density and with 
power law slope — p (not to be confused with pressure p). The 
fraction of available electrons that is accelerated is given by 
A further fraction eg of the local internal energy density 
resides in the shock-generated magnetic field. The effect of 
electron cooling is included using a global estimate for the 
electron cooling time t c , by equating it to the lab frame time 
since the explosion. The spectral shape of the absorption coef- 
ficient a v for s.s.a. is also given by sharply connected power 
laws and the effect of electron cooling on ot v is ignored (in 
any case, the error thus introduced is neglible compared to 
the small error from using a global rather than a local esti- 
mate for electro n cooling. Global and lo cal electron cooling 
are compared in [Van Eerten et al.ll2010al) . Mathematical ex- 
pression s for the emission and absorption coefficients can be 
found in lVan Eerten et al1l2012l 

After all rays have emerged from the blast wave and the ob- 
served flux is calculated by integrating over the rays, the ob- 
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served spectrum will again consist of a set of power laws, now 
smoothly connected due to the different break frequencies at 
different contributing parts of the fluid. S.s.a. manifests itself 
as an additional break in the spectrum, occurring typically at 
radio wavelengths. An example spectrum, calculated at early 
time for Qq = 0.05 rad and ISM environment is provided by 
Fig. 12 A comparison bet ween a light curve from the boosted 
frame and from BOXFIT (IVan Eerten et al.ll20l2l) is given by 
the bottom plot of Fig. Q] Although Fig. [2] includes s.s.a., 
we will focus in this work on optical and X-ray frequencies 
(where the jet break is typically observed) and postpone a de- 
tailed treatment of the s.s.a. break to future work. 

3.1. Resolution 

All light curves were calculated for p = 2.5, which is suf- 
ficient to derive light curves for arbitrary p > 2 value, as ex- 
plained below in section |4] Lacking an upper cut-off to the 
accelerated particle energy distribution, our radiation model is 
invalid for p <= 2 because then the integral for the total accel- 
erated particle energy diverges. The following settings were 
used to compute numerically converged synthetic light curves 
from the boosted frame simulations. 3000 simulation snap- 
shots were probed for computing on-axis light curves, 300 
snapshots were probed for computing off-axis light curves. A 
matrix of rays was used consisting of 1500 rays logarithmi- 
cally spaced in the radial direction and 100 evenly spaced in 
the angular direction (or 1 in the angular direction, for on- 
axis observations). These directions refer to coordinates on 
the plane perpendicular to the observer (in the lab frame, see 
appendix[A}, and the inner and outer boundaries of this plane 
are 10 12 and 10 18 cm respectively. Each light curve has 150 
data points between observer times of 10~ 4 days and 10 2 days, 
although only 85 (60) data points between 0.01 (0.1) and 26 
days were used for analysis in order to ensure complete cover- 
age of the observer times by the emission times. Light curves 
were calculated for half-opening angles 6q = 0.05, 0.1 and 0.2 
rad and observer angles 9 / n that were a fraction 0, 0.2, 0.4, 
0.6, 0.8 or 1 of 6>o. 

4. PEAK FLUX AND BREAK FREQUENCIES 
4. 1 . Theory 

In section |2~21 we demonstrated the scale-invariance of the 
jet dynamics between different jet energies and circumburst 
densities. In [Van Eerten & MacFadven (12012a!) we demon- 
strated that similar scalings apply in the asymptotic regimes 
of the observed spectrum. Although additional constants are 
introduced when calculating synchrotron radiation, such as 
the electron mass m e , these can be identified and isolated in 
the flux equations in a given spectral regime and the resulting 
scale invariance is again a result of dimensional analysis. 

Invariance of the fluxes in each asymptotic spectral regime 
is equivalent to scale invariance of the critical quantities that 
determine the shape of the spectrum: peak flux F pea k, syn- 
chrotron break frequency u m , cooling break frequency v c and 
s.s.a. break frequency v a . The shapes of the spectral tran- 
sitions are not scale-invariant, but can be m odeled using a 
smoo t h connection between power laws (see iGranot & Saril 
12001 IVanEerten & Wiiersl 120091: ILeventis et al.ll2012l) . The 
flux for a given observation can be calculated from the charac- 
teristic scale-invariant evolution of the critical quantities plus 
a description of the spectral transitions (which can also be a 
simple sharp power law approximation), and from three sets 
of parameters defining the observation and model: the ob- 



10 53 erg J 

T = (\/K) X ^- k) t ohs /(\+Z) 



(\ + Z) p— 1 n 1/2 i 3<2-t) 3 

d\ % 3/7-1 



(l + z) p-l o l/2 t l \l/2 £ / a a l n\ 
Fpeak= 2 ? _-l € e e B ?JV kA ]peak(T\ , obs , k = 0) 
"28 

v m = (l+zT l (jZj) e 2 e e B /2 ^ 2 K X l / 2 f m (T-6 A,bs,k = Q) 
v c = (l + zr l ey B y2 ^K- 2 ' 3 \- 5 ^ c (T-6 ,9 obs ,k = 0) 



Fpeak = —n— ir — Te^eB 2 ^l/K X } ^ 2 fp ea k(T;9o,0 h s ,k = 2) 

v m = {l+zT l e 2 e B /2 ^N^ 1 ^ /2 Ur\doA>b s ,k = 2) 

u c = (l + zT\° e e B 3/2 ^ NK 1 y 5/2 fc(T-,Oo,e o bs,k = 2) 



TABLE 1 

Concise equations for the characteristic quantities. Top 
panel shows energy and density scale factors k and a and 
scaled time t. the following panels show equations for 
general k, k = (ism) and k = 2 (stellar wind) from top to 

BOTTOM. 

server parameters z (redshift), cIl (luminosity distance), o b s 
(observer angle), t a \, s (time) and v (frequency); the explosion 
parameters k, A, Ei SO , 9q; the radiation parameters p, eg, e e , 
£ai. The dependency of the flux on d.L and z is straightforward, 
with F v oc dj 2 and flux, frequency and time depending on z in 
the standard manner. Different values for 9„b s , &o and k lead 
to different evolution of the characteristic quantities v a , v m , v c 
and F pea k. Scale invariance takes care of E iso and A, while the 
dependency of the characteristic quantities on eg, e e , £jv re- 
mains unchanged throughout the evolution of the decollimat- 
ing blast wave and can be determined analytically for general 
p and k. The dependency of the characteristic quantities on p 
is constant in time and known analytically, meaning that once 
the evolution for a given p value is known, their evolution for 
any p value can be trivially obtained. Different spectral order- 
ings of the break frequencies v a , v m , v c lead to different time 
evolutions. 

Table Q] summarizes these properties of the light curves. 
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Here the general k case is shown as well as the ISM and 
stellar wind cases separately. o?28 is the luminosity distance 
di in units of 10 28 cm. The functions f pea k(T',9o,9 b s ,k), 
fm(T",9o,d b s ,k) and f c (r; 9o,6 b s ,k) denote the scale-invariant 
time evolution of the characteristic quantities as they are de- 
termined numerically from analyzing light curves computed 
from the boosted frame simulations for each spectral regime. 
These functions can be scaled from their baseline values to 
arbitrary explosion energy and circumburst density by plug- 
ging in the scaled values for k, A and r from the top section 
of the table into the equations in the lower sections of the ta- 
ble. Since their dependency on the radiation parameters eg, 
e e , £n and p is known and constant in time, these terms have 
been made explicit in the table. Left implicit is the fact that 
different spectral orderings lead to different evolution curves, 
which will affect the characteristic scale-invariant functions f 
but not their pre-factors. In the remainder of this paper we 
will discuss the slow cooling case (which is the one usually 
observed in practice), where v m < v c . The observed flux fol- 
lows from the characteristic quantities according to 



around 



tj w (0.6±0.1)(1+z)(k/A) i/3 (6» /0.1) 8/3 days, 



(5) 



Fd - Fpeak 



Fg — Fpeafc ( 



1/3 



d-p)/2 



V < v m < v c , 



V m < V < V c , 



F H = F peak [—J I — J > v m < v c < v, (4) 



where the l abels D, G, H have b een chosen to match the no- 
tation from lGranot & Saril (120021) . Note that, even though the 
time evolution of the characteristic quantities does not depend 
on p, equations|4]implies that the time evolution of Fc and Fh 
does. 



4.2. Numerical results 

In Fig. [3] we plot the time evolution of the characteristic 
quantities for the three jet opening angles 9q = 0.05, 0.1, 0.2 
rad. and for observer angles 6 t, s = 0, 0.6 x 9q, 9q. Fig. |4]shows 
the evolution of the spectral slope for these same angles. The 
time evolutions were calculated from three light curves per 
9o, dobs combination: one for each asymptotic spectral regime 
separated by v m and v c in the slow cooling case. For these 
light curves we used cb = 10~ 5 , e e = 10~ 5 , = 1, p = 2.5 and 
frequencies 10~ 20 , 10 10 , 10 40 Hz. These values (especially 
the frequencies) were not physically motivated but rather cho- 
sen such that they ensure all light curves were calculated well 
into the asymptotic limits of the spectral regimes and not im- 
pacted by the smoothness of the spectral transitions between 
regimes. Moving the outer frequencies closer in but still in 
their asymptotic regions throughout the evolution of the emis- 
sion for eg = 10~ 5 , e e = 10~ 5 , = T P = 2.5, e.g. to 10~ 5 and 
10 25 was found to have no impact on the result. The char- 
acteristic quantities were subsequently obtained by inverting 
equations!?] Since the plots show f pea k, f m and f c , the curves 
are independent of e e , eg and p. 

Figs. |3] and |4] reveal that the time evolution for the charac- 
teristic quantities strongly depends on both jet and observer 
angle. A difference between pre- and post-break values sets 
in immediately following the jet break time tj, which for the 
current ISM simulations and on-axis observers is found to lie 



consistent with t he earli er reported numerical results from 
IVan Eerten et al.l (l2010al) . The jet break time will be dis- 
cussed in more detail in section[5]below. In Fig. @] both pre- 
and post-break theoretically expected slopes are also plotted. 
The pre-break slopes match the theoretical predictions well, 
but the post-break s lopes differ su bstantially from theoreti- 
cal predictions (see lSari et al.lll999tlRho ads 1999). This is a 
consequence of the fact that the spreading of blast waves in 
simulations is not a runaway process even for 9q = 0.05 rad. 
jVan Eerten & M acFadvenl l2012bt iMacFadven & van Eertenl 
120121) . 

It would be a strong indication of self-similarity between 
jet opening angles and of great practical significance if the 
evolution functions were to scale in a straightforward manner 
between opening angles. Any such scaling should incorporate 
the 9q dependency of jet break time tj. When we take that as 

our starting point, scale time according to t' = t(9' /9o) s ' 3 and 
the characteristic functions according to f(f') = (9' /9o)~ a f(t), 
where a is the power law time-dependence in the pre-break 
BM regime, this yields evolution curves for f pea ^ that numer- 
ically match quite well initially between different jet open- 
ing angles, even for off-axis observer angles, as illustrated in 
Fig. [5] These #o-scalings however are not exact and a similar 
mapping for f„, or f c fails to produce much numerical overlap. 
This can be seen from the power law slope plots in Fig. [4] 
since the scalings represent horizontal shifts of the curves in 
these plots. Essentially, this lack of straightforward scalabil- 
ity reflects the fact that the post-break behavior is determined 
by more characteristic timescales than tj alone, such as the 
transition time to non-relativistic flow and the transition time 
to quasi-spherical flow, and that these timescales will impact 
the trans-relativistic stage of fluid flow that generates the ob- 
served post-break light curves. 

4.3. The cooling break 

Figs. [3] and |4] illustrate that for a given characteristic func- 
tion different extremal values are reached for different ob- 
server and jet angles. The consequence of this on the jet 
break as measured from observations will be discussed in 
section [5] below, and we limit ourselves here to highlight- 
ing the behavior of the cooling b reak. In an earlier work 
dVan Eerten & M acFadvenl l2012al) we showed how simula- 
tions in a fixed frame (and thus of lower resolution) indicated 
a steepening (for 9„b s = 0) of the temporal evolution of the 
cooling break immediately following the jet break. However, 
the current boosted frame simulations reveal no post-break v c 
steepening towards stronger decay. Specifically, the v c light 
curves for 9p = 0.2 rad., th e same a ngle as plotted in Fig. 3 of 
IVan Eerten & MacFadven (2012a), show only a turnover to- 
wards positive temporal slope following the jet break, as can 
be seen in the bottom panel of Fig. |4]of the current paper. At 
the same time, the on-axis curve for 9 b s = 0.05 rad. shown in 
the same figure, does show a (slight) post-break steepening of 
the temporal power law slope of v c . 

What this indicates is that the earlier reported steepen- 
ing for 9q = 0.2 rad. and the current smaller steepening 
for 9o = 0.05 rad. are numerical in origin, and sensitive to 
the initial conditions of the blast wave. Above the cool- 
ing break, the observed flux is dominated by emission from 
the edges of the jet (i.e. the observed image is 'limb- 
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FIG. 3. — Time evolution of characteristic quantities, for different jet open- 
ing angles and observer angles. Top to bottom: f pea t, f,„, f c . They describe 
the time evolution of Fp^, v m , u c respectively according to the equations in 
table [T] The legend in the top plot applies to all plots. 



brightened'), and therefore the cooling break v c is the most 
sensitive to deviations from purely radial flow at the edges 
of an initially conically truncated spherical BM outflow. The 
smaller the jet opening angle, the larger even small resolution- 
induced deviations are relative to 6*o- The dynamics of narrow 
and wide jets will be discussed separately in more detail in 



FIG. 4. — Time evolution of power law slopes of the characteristic quan- 
tities, for different jet opening angles and observer angles. Top to bottom: 
slopes for f pea k, f m , f c . The legend in the top plot of Fig. reapplies to all these 
plots as well. The constant grey lines indicate the expected slopes for light 
curves for an on-axis observer from the pre-break BM solution and assuming 
a fast spreading jet post-break. 
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10" 1 10° 10 1 

r (days) 

FIG. 5. — Scaled evolution of the peak flux function f pea t, where the curves 
for 8q = 0.05 rad. and 8q = 0.2 rad. have been scaled in time towards the 
6»o = 0.1 rad. result using t' = (8 /0.lf/ 3 . As in Figs. |5]and|4] solid lines 
refer to 8 \, s = 0, dotted lines to 8„i, s = 0.68q and dashed lines to 8 t, s = 8q. 



well to the left of the left vertical line at 10~ 2 days. In view of 
this resolution issue, when analysing light curves for do = 0.2 
rad we will start from 0.1 days (rather than 0.01 days), and 
the characteristic evolution curves for this initial jet opening 
angle in Figs|3]and|4]have been truncated at this value of r. 

An additional conclusion that can be drawn from the severe 
resolution dependence of the off-axis observed v c evolution 
for #o = 0.2 rad., even well into times that are easily observ- 
able by instruments such as Swift, is that if small numerical 
resolution-induced deviations from BM-type flow will have a 
large effect on v c , the same will hold for minor physical devia- 
tions. This renders relevant the question to what extent devia- 
tions from the expected BM-based time evolution of v c can be 
driven by the dynamics of the outflow. On the other hand, al- 
though an actual measurement of the evolution of v c has been 
performed by Filgas et al.J (1201 11) . the temporal slope of -1.2 
that these authors find is steeper than the high-resolution sim- 
ulation v c slope in Fig. [4] at any time, and their thesis is that 
the steep decline in GRB 091 127 can be attributed to changes 
in the radiative process (via a time dependency of eg) rather 
than outflow dynamics. 
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FIG. 6. — A comparison of the time evolution of four computations of v c 
for 8q = 0.2 rad. and 8 ], s — 8q (i.e. an on-edge observer). The vertical lines 
indicate the timespans used for analysis. 



Mac Fadven & van Eerteiil (120121) . 

The effect of early time flow at the jet edges on the light 
curve naturally becomes more severe the closer the observer 
angle moves towards the edge of the jet. In Fig. [6] we show 
that even the pre-break behavior for v c from 2D simulations 
differs strongly from that from analytically calculated coni- 
cal outflow. The black solid line and blue dashed line show 
v c results for simulations starting at 70 = 100 and 70 = 50 re- 
spectively, the red dash-dotted line shows the evolution of v c 
based on conical outflow following the BM solution, while the 
green dashed line shows a four times lower resolution simu- 
lation. Around the leftmost vertical line, the v c curves for 
both normal resolution simulations have merged, while both 
simulation curves still differ strongly from the BM solution. 
It follows that the difference between 2D simulated and ra- 
dial analytical flow can not be attributed to a lack of early 
time coverage of the observed signal by an incomplete range 
of emission times. Nor can this difference be attributed to 
the difference in starting times (and hence the extent to which 
70 <C l/#o- Both effects are clearly visible in Fig [6] and lie 



5. LIGHT CURVES AND JET BREAKS 

Once the time evolution of the characteristic functions f pea k, 
f c and f c is known, they can be used to quickly calculate light 
curves for arbitrary p. In order to study the shape of the jet 
break we have done so for p values of 2.01 and 2.1, 2.2, . . ., 
3.0 and the spectral regimes v a < v m < v < v c (typically appli- 
cable to optical data) and v a < v m < v c < v ("X-rays"). Vari- 
ous functions have been used in the literature to fit jet breaks 
in opti cal and X-ray light curve data, such as s harp power laws 
(e.g. iRacusin et al J I2009t lEvans et al.l 12 009). smoothly con- 
nected power laws (e.g. iBeuermann et alJI 1999b or power law 
transit ions where the turn over includes an exponential term 
(e.g. lHarrison et alj [l999). A limitation common to all these 
fit functions is the assumption of a single power law regime 
after the jet break. Although simulation based light curves 
show that in reality this should not be expected to be the case 
(as can be seen from the post-break evolution of the peak flux 
and break frequencies, Figs. [3]and|4|i, it therefore makes prac- 
tical sense to explore the implications of our simulation re- 
sults for the interpretation and applicability of broken power 
law fit functions. 

Tables [2] and [3] show the results of the analysis of light 
curves with different p, 8 i, s and 6$ values using different 
power law descriptions. Each data point of the light curves, 
consisting of 85 data points per curve between observer times 
10~ 2 days and 26 days for Oq = 0.05 rad. and 80 = 0.1 rad. and 
60 per curve between observer times 10 _1 days and 26 days 
for 6*o = 0.2 rad., was given an error of ten percent and three 
different jet break functions were fitted using a least squares 
algorithm. A baseline frequency v = 4.56 x 10 14 Hz (R-band) 
was used for tableland a baseline frequency v = 5 x 10 17 Hz 
(2.07 KeV) for table[3] The fit function for a sharp power law, 
labeled "PL" in the table and below, is given by 



F(r): 



C{T/T b ) a \ T<T h , 

C{T/nT\ T >n 



(6) 



The fit function for a smooth power law transition, equivalent 
to that used by Be uermann et al.l [1999 and labeled "sB", is 
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6q (radians) o i, s 


fit function 


a 


Oil 


< T A > 


10 log C 


<7 






X , reduced^ 


0.05 


PL 


0.76-0.73/) 


0.20-1.24/) 


0.10 


1.16 + 2.27/) 






1 


0.20 




sB 


0.75-0.71/) 


0.19-1.24/) 


0.10 


1.13 + 2.30/) 


6.06 


-0.91/) 


0.68 


0.14 




sH 


1.01-0.74/; 


0.18-1.23/) 


0.08 


1.36 + 2.29/) 






1.00 


0.20 


0.26*0 


PL 


0.69-0.74/) 


0.19-1.24/) 


0.12 


1.13 + 2.23/) 






1 


0.54 




sB 


0.71-0.69/> 


0.20-1.25/) 


0.10 


1.12 + 2.31/) 


3.65 


-0.70/) 


0.34 


0.18 




sH 


1.02-0.78/) 


0.19-1.24/) 


0.09 


1.45 + 2.21/) 






0.28 


0.15 


0.46*0 


PL 


0.67-0.82p 


0.23-1.28/) 


0.17 


1.33 + 1.98/) 






1 


0.99 




sB 


0.72-0.77/) 


0.23-1.29/) 


0.15 


1.32 + 2.08/) 


2.92 


-0.51/) 


0.58 


0.57 




sH 


0.94-0.85/) 


0.22-1.28/) 


0.14 


1.65 + 1.98/) 






0.46 


0.46 


0.66*0 


PL 


0.56-0.80/) 


0.26-1.32/) 


0.25 


1.31 + 1.85/) 






1 


0.74 




sB 


0.58-0.79/) 


0.26-1.32/) 


0.24 


1.31 + 1.87/) 


5.63 


-0.92/) 


0.85 


0.64 




sH 


0.78-0.83/) 


0.24-1.32/) 


0.21 


1.62+1.82/) 






0.87 


0.65 


0.86*0 


PL 


0.56-0.77/) 


0.29-1.37/) 


0.33 


1.23 + 1.79/) 






1 


1.07 




sB 


0.56-0.75/) 


0.29-1.37/) 


0.32 


1.19+1.83/) 


5.72 


-1.09/) 


0.86 


0.93 




sH 


0.74-0.79/) 


0.27-1.37/) 


0.29 


1.48 + 1.77/) 






0.84 


089 


e 


PL 


0.66-0.77/) 


0.31-1.41/) 


0.43 


1.13 + 1.75/) 






1 


1.48 




sB 


0.68-0.76/) 


0.33-1.42/) 


0.41 


1.18 + 1.76/) 


4.26 


-0.74/) 


0.82 


1.23 




sH 


0.82-0.79/) 


0.31-1.42/) 


0.38 


1.41 + 1.72/) 






0.75 


1.11 


0.1 


PL 


0.74-0.75/) 


0.26-1.29/) 


0.59 


1.69+1.70/) 






1 


0.23 




sB 


0.75-0.75/) 


0.25-1.29/) 


0.58 


1.70+1.71/) 


5.35 


-0.77/) 


0.61 


0.14 




sH 


0.85-0.76/) 


0.21-1.28/) 


0.54 


1.85 + 1.70/) 






1.11 


0.24 


O.20o 


PL 


0.73-0.76/) 


0.26-1.29/) 


0.64 


1.71 + 1.66/) 






1 


0.66 




sB 


0.74-0.74/) 


0.26-1.31/) 


0.63 


1.72+1.70/) 


2.96 


-0.52/) 


0.17 


0.12 




sH 


0.85-0.77/) 


0.22-1.29/) 


0.59 


1.90+1.65/) 






0.10 


0.06 


0.46*0 


PL 


0.73-0.80/) 


0.32-1.34/) 


0.91 


1.93 + 1.44/) 






1 


2.06 




sB 


0.75-0.74/) 


0.30-1.40/) 


0.86 


1.91 + 1.56/) 


1.49 


-0.27/) 


0.20 


0.42 




sH 


0.85-0.81/) 


0.27-1.34/) 


0.82 


2.11 + 1.44/) 






0.29 


0.62 


0.66*0 


PL 


0.67-0.82/) 


0.44-1.44/) 


1.45 


1.95 + 1.25/) 






1 


1.90 




sB 


0.74-0.81/) 


0.23-1.43/) 


1.52 


2.10+1.25/) 


1.22 


-0.12/) 


0.46 


0.89 




sH 


0.75-0.82/) 


0.32-1.42/) 


1.36 


2.17 + 1.22/) 






0.44 


0.84 


0.86»o 


PL 


0.59-0.79/) 


0.43-1.49/) 


2.19 


1.78 + 1.17/) 






1 


0.59 




sB 


0.59-0.78/) 


0.42-1.54/) 


2.29 


1.87 + 1.15/) 


3.20 


-0.59/) 


0.28 


0.17 




sH 


0.66-0.80/) 


0.28-1.48/) 


2.19 


1.97 + 1.13/) 






0.18 


0.10 


do 


PL 


0.71-0.79/) 


0.51-1.56/) 


2.67 


1.85 + 1.09/) 






1 


1.30 




sB 


0.71-0.77/) 


0.46-1.66/) 


3.07 


1.88 + 1.08/) 


1.90 


-0.35/) 


0.27 


0.36 




sH 


0.76-0.80/) 


0.34-1.55/) 


2.74 


1.95 + 1.07/) 






0.24 


0.33 


0.2 


PL 


0.79-0.80/) 


0.51-1.44/) 


3.73 


2.40+1.03/) 






1 


0.16 




sB 


0.78-0.79/) 


0.45-1.44/) 


3.77 


2.37 + 1.04/) 


5.89 


-0.96/) 


0.36 


0.06 




sH 


0.88-0.80/) 


0.25-1.39/) 


3.64 


2.47 + 1.04/) 






1.57 


0.23 


0.26*0 


PL 


0.79-0.81/) 


0.53-1.40/) 


3.75 


2.42+1.00/) 






1 


0.59 




sB 


0.78-0.78/) 


0.44-1.48/) 


4.25 


2.41 + 1.01/) 


2.69 


-0.50/) 


0.05 


0.03 




sH 


0.88-0.82/) 


0.23-1.36/) 


3.87 


2.51 + 1.00/) 






0.03 


0.02 


0.46»o 


PL 


0.77-0.85/) 


0.57-1.36/) 


4.54 


2.51+0.87/) 






1 


1.71 




sB 


0.83-0.78/) 


-0.43 -2.22p 


21.03 


3.01+0.49/) 


0.34 


-0.05/) 


0.06 


0.11 




sH 


0.88-0.86/) 


0.32-1.36/) 


4.86 


2.69 + 0.82/) 






0.25 


0.44 


0.66*0 


PL 


0.70-0.86/) 


0.80-1.45/) 


7.38 


2.69 + 0.61/) 






1 


1.07 




sB 


0.74 - 0.85/) 


0.70-2.55/) 


24.27 


2.93 + 0.28/) 


0.60 


-0.10/) 


0.28 


0.31 




sH 


0.76-0.87/) 


0.18-1.37/) 


8.67 


2.65 + 0.62/. 






0.36 


0.39 


0.86*0 


PL 


0.62-0.83/) 


0.64-1.26/) 


9.42 


2.41+0.62/) 






1 


0.26 




sB 


0.60-0.81/) 


0.15-2.14/) 


24.47 


2.67 + 0.30/) 


1.25 


-0.26/) 


0.05 


0.01 




sH 


0.61-0.81/) 


-l.38-0.86/) 


14.05 


2.14 + 0.65/) 






0.05 


0.01 


Oo 


PL 


0.73-0.81/) 


0.63-1.03/) 


6.96 


2.08 + 0.82/) 






1 


0.25 




sB 


0.71-0.79/) 


0.21-1.50/) 


24.46 


2.58 + 0.36/) 


1.44 


-0.29/) 


0.14 


0.03 




sH 


0.72-0.80/) 


-0.80-0.77/) 


12.83 


1.90 + 0.78/) 






0.13 


0.03 



TABLE 2 

FIT RESULTS FOR U obs < V c . PL, sB AND sH REFER TO DIFFERENT FIT FUNCTIONS DESCRIBED IN THE TEXT AND Ct , Cti , T b , C AND a TO THEIR 
PARAMETERS. ^ PLEASE NOTE THAT THE SCALE OF THE reduced \ 2 RESULT IS ARBITRARY BECAUSE IT DEPENDS ON THE NUMBER (85 OR 60) AND 
SPACING (LOGARITHMIC IN TIME) OF SYNTHETIC LIGHT CURVE DATAPOINTS (CHOSEN PURELY TO PROPERLY RESOLVE THE LIGHT CURVE EVOLUTION) 
AS WELL AS AN ARBITRARY ERROR OF 10 % ON EACH DATAPOINT, WHICH HAVE NO PHYSICAL MEANING. 
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9q (radians) 


dobs 


fit function 


ao 


ai 


< 7J > 


10 log C 




x 2 /x 2 p, 


X 2 , reduced 


0.05 





PL 


0.50-0.72p 


0.20-1.23/7 


0.08 


-0.24 + 0.83/7 




1 


1.92 






sB 


0.50-0.71/7 


0.20-1.23/7 


0.07 


-0.25 + 0.83/7 


23.81-4.90/7 


1.00 


1.94 






sH 


0.93-0.76/7 


0.20-1.23/7 


0.06 


0.13 + 0.81/7 




1.09 


2.10 




O.20 o 


PL 


0.50-0.74/7 


0.21-1.23/7 


0.08 


-0.18 + 0.77/7 




1 


1.86 






sB 


0.49-0.71/7 


0.22-1.24/7 


0.08 


-0.22 + 0.81/7 


10.25-2.27/7 


0.97 


1.82 






sH 


0.95-0.80/7 


0.21-1.24/7 


0.06 


0.22 + 0.75/7 




0.98 


1.83 




0.46*0 


PL 


0.47-0.82/7 


0.26-1.27/7 


0.12 


0.03 + 0.51/7 




1 


2.29 






sB 


0.52-0.79/7 


0.27-1.27/7 


0.11 


0.03 + 0.58/7 


5.26-0.99/7 


0.92 


2.14 






sH 


0.89-0.88/7 


0.27-1.27/7 


0.09 


0.51+0.49/7 




0.88 


2.00 




0.66*0 


PL 


0.33-0.80/7 


0.33-1.31/7 


0.20 


-0.06 + 0.35/7 




1 


2.17 






sB 


0.33-0.80/7 


0.33-1.31/7 


0.19 


-0.06 + 0.37/7 


13.33-2.63/7 


0.99 


2.17 






sH 


0.66-0.85/7 


0.33-1.32/7 


0.14 


0.46 + 0.30/7 




1.06 


2.31 




0.86*0 


PL 


0.32-0.77/7 


0.39-1.36/7 


0.27 


-0.19 + 0.30/7 




1 


2.37 






sB 


0.31-0.75/7 


0.39-1.36/7 


0.26 


-0.23 + 0.33/7 


11.95-2.67/7 


0.98 


2.34 






sH 


0.60-0.81/7 


0.39-1.37/7 


0.21 


0.33 + 0.22/7 




1.03 


2.44 




Oo 


PL 


0.45-0.78/7 


0.46-1.42/7 


0.35 


-0.27 + 0.23/7 




1 


2.55 






sB 


0.43-0.76/7 


0.46-1.42/7 


0.33 


-0.29 + 0.27/7 


7.98-1.67/7 


0.96 


2.47 






sH 


0.66-0.81/7 


0.45-1.42/7 


0.29 


0.13 + 0.18/7 




0.97 


2.45 


0.1 





PL 


0.50-0.75/7 


0.27-1.29/7 


0.52 


0.11+0.22/7 




1 


0.60 






sB 


0.50-0.74/7 


0.28-1.29/7 


0.51 


0.10 + 0.23/7 


12.62-2.52/7 


0.95 


0.58 






sH 


0.66-0.77/7 


0.25-1.28/7 


0.44 


0.39 + 0.19/7 




1.40 


0.84 




0.26*0 


PL 


0.49-0.75/7 


0.27-1.29/7 


0.55 


0.09 + 0.20/7 




1 


0.65 






sB 


0.49-0.74/7 


0.30-1.31/7 


0.53 


0.13 + 0.21/7 


5.96-1.26/7 


0.67 


0.44 






sH 


0.65-0.78/7 


0.26-1.29/7 


0.48 


0.43 + 0.14/7 




0.69 


0.44 




0.46*0 


PL 


0.49-0.80/7 


0.34-1.32/7 


0.76 


0.30-0.03/7 




1 


1.64 






sB 


0.49-0.74/7 


0.41-1.39/7 


0.68 


0.30+0.07/7 


2.53-0.53/7 


0.38 


0.63 






sH 


0.66-0.82/7 


0.33-1.33/7 


0.65 


0.67 - 0.08/7 




0.33 


0.55 




0.66*0 


PL 


0.43-0.81/7 


0.48-1.42/7 


1.26 


0.42-0.28/7 




1 


1.50 






sB 


0.48-0.81/7 


0.44-1.44/7 


1.25 


0.49-0.27/7 


2.22-0.32/7 


0.63 


0.96 






sH 


0.56-0.84/7 


0.43-1.42/7 


1.14 


0.72-0.33/7 




0.50 


0.76 




0.86*0 


PL 


0.35-0.79/7 


0.54-1.49/7 


2.00 


0.17-0.36/7 




1 


0.40 






sB 


0.34-0.78/7 


0.59-1.53/7 


2.04 


0.19-0.36/7 


6.73-1.53/7 


0.56 


0.22 






sH 


0.44-0.81/7 


0.42-1.48/7 


1.96 


0.36-0.40/7 




0.81 


0.30 




8o 


PL 


0.48-0.79/7 


0.71-1.58/7 


2.52 


0.23-0.45/7 




1 


0.82 






sB 


0.46-0.77/7 


0.74-1.66/7 


2.69 


0.16-0.42/7 


3.48-0.76/7 


0.39 


0.32 






sH 


0.54-0.80/7 


0.49-1.54/7 


2.50 


0.28-0.44/7 




0.28 


0.23 


0.2 





PL 


0.56-0.80/7 


0.52-1.45/7 


3.66 


0.62-0.49/7 




1 


0.14 






sB 


0.55-0.79/7 


0.49-1.44/7 


3.65 


0.59-0.46/7 


9.44-1.85/7 


0.61 


0.09 






sH 


0.67-0.80/7 


0.28-1.40/7 


3.48 


0.74-0.48/7 




2.72 


0.36 




0.26»o 


PL 


0.57-0.81/7 


0.53-1.41/7 


3.65 


0.66-0.51/7 




1 


0.40 






sB 


0.54-0.78/7 


0.54-1.48/7 


3.96 


0.63-0.50/7 


4.07-0.85/7 


0.10 


0.04 






sH 


0.67-0.82/7 


0.26-1.37/7 


3.67 


0.77-0.52/7 




0.15 


0.05 




0.46*0 


PL 


0.57-0.85/7 


0.66-1.39/7 


4.28 


0.90-0.70/7 




1 


1.28 






sB 


0.58-0.78/7 


1.21-2.55/7 


18.24 


1.26-1.07/7 


0.55-0.11/7 


0.06 


0.08 






sH 


0.68-0.86/7 


0.33-1.35/7 


4.55 


0.97-0.70/7 




0.19 


0.26 




0.66*0 


PL 


0.46-0.86/7 


0.62-1.36/7 


7.33 


0.75-0.85/7 




1 


0.72 






sB 


0.51-0.85/7 


1.03-2.48/7 


24.07 


0.89-1.20/7 


0.74-0.13/7 


0.29 


0.21 






sH 


0.54-0.87/7 


0.12-1.37/7 


9.25 


0.75-0.89/7 




0.33 


0.24 




0.86*0 


PL 


0.40-0.82/7 


0.45-1.21/7 


9.71 


0.43-0.85/7 




1 


0.15 






sB 


0.39-0.81/7 


0.42-2.14/7 


24.29 


0.71-1.21/7 


1.63-0.36/7 


0.09 


0.01 






sH 


0.39-0.81/7 


-2.03-0.72/7 


16.11 


0.10-0.84/7 




0.09 


0.01 




8o 


PL 


0.50-0.81/7 


0.42-1.01/7 


7.82 


0.11-0.67/7 




1 


0.13 






sB 


0.48-0.80/7 


0.12-1.48/7 


24.44 


0.62-1.16/7 


2.03-0.45/7 


0.08 


0.01 






sH 


0.49-0.80/7 


-1.70-0.58/7 


16.48 


-0.14-0.74/7 




0.06 


0.01 



TABLE 3 

Same as table[2] now for v obs > v c . 
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given by 



F(t) = C 



-I/O 



(7) 



The alternative smooth power law transition fit function, la- 
beled "sH" is the same as the one used by lHarrison et al.l 
( 119991) and given by 



F(T) = C{l-CXp[-(T/T h ) a "- a ']}(T/T b r' 



(8) 



In this equation the pre-break power law slope is retrieved 
from the Taylor series of the exponential term. 

The different fit variable results are represented in the tables 
as follows. Since the fit results confirm that the slopes ao and 
a\ linearly depend on p (as shown in Figs. [7J and [8] for fits 
using sharp power laws), the entries contain this linear depen- 
dence as determined from the full range of p fits rather than 
a values for each individual p. The logarithm of the numeri- 
cal scale factor C and the sharpness of the smooth power law 
fit also depend linearly on p, and are presented in the same 
fashion. The break time r depends only weakly on p and is 
represented by its average value < r >, weighing equally all 
individual p value fits. We also give the reduced \ 2 value of 
each fit, again averaged over the different p value fits, as well 
as the ratio of the unreduced \ 2 of each fit to that of a sharp 
power law fit. For the latter, these ratios were calculated be- 
fore the average was taken. We emphasize that by themselves, 
the reduced \ 2 results are to some degree arbitrary, since they 
depend on arbitrary quantities like the number of datapoints 
in a synthetic light curve, the spacing of these datapoints and 
an artificial ten percent error on each datapoint, and that they 
should be interpreted only in a relative sense. 

Using the prescriptions from table Q] and equations [4] the 
flux for any combination of parameter values can be repro- 
duced from the fit results in tables [2] and [3] For v < v c , we 
get: 



(l +z )(3-p)/2 p _ l / p _ 2 



d 2 

"28 



/'-I 



K l A (p+D/4 



3/7-1 \p-l 
4.56 x 10 14 Hz 



(1-/0/2 



F(t) mJy, (9) 



for the ISM case and F referring to fit results from table[2] We 
have now added a '©' to the frequency to emphasize that this 
frequency is expressed in the observer frame, like the char- 
acteristic frequencies in table Q] and the frequencies in Eq. 
|H and related to the frequency v in the burster frame via 
v ® = vj(\ +z). Note that r is still expressed in the burster 
frame, in order to keep the characteristic functions (and hence 
the power law fit results) redshift-independent. For v > v c we 
have: 



F H = 



(l +z )(2-rt/2 p _ l / p _ 2 



d 2 

"28 



3p-l \p-l 



p-l (p-2)/4 f .2-p 
SA' 



K 2/3 A (3p-2)/12 



-p/2 



F(t) mJy, (10) 



s 5 x 10 17 Hz / 

for the ISM case and F referring to fit results from table [3] 

5.1. Implications for the light curve slope 

The predicted on-axis pre-break slopes are 3(1 — p)/4 for 
v < v c and (2-3/?)/4 for v > v c and the tables show that 
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FIG. 7. — Pre-break temporal index ao (top plot) and post-break temporal 
index ot\ (bottom plot) for v < v c and for on-axis and on-edge observations 
of the three jet angles 6q = 0.05, 0.1, 0.2 rad., according to sharp power law 
fits to synthetic light curves. The grey bands indicate the region within ten 
percent of the theoretically expected values, 3(1 — p)/4 and — p for pre- and 
post-break respectively. 



these values are well reproduced by straight power law fits 
for 9q = 0.05 rad. and 8q = 0.1 rad. and reasonably well for 
6*o = 0.2 rad. This can also be seen from the top plots in Figs. [7] 
and[8] that show ao for each p value and on-axis and on-edge 
observers. 

The post-break slopes, on the other hand, are not con- 
sistent with the theoretically expected temporal index —p 
Sari et all 1 1999) nor with very smooth gra d ual tr ansitions 



Kumar & Panaitescu 200Qt iWei & Lul l2000l l2002h . as can 
be seen from the tables and the bottom plots in Figs. [7] and 
[8] and are steeper to the extent that they fall well outside 
even a ten percent margin of the theoretical value. We find 
a steepening of about -(0.5+0.5/?) below the cooling break 
and about -(0.25 + 0.5/?) above the cooling break, both lead- 
ing to post-break slopes of roughly 0.25 — 1.3/?, although dif- 
ferent observer angles, jet opening angles and heuristic de- 
scriptions of the break introduce a wide range of temporal 
indices. This confirms earlier numerical work, and was first 
shown fro m high-resolution simulation s and for on-axis ob- 
servers by Zhang & MacFadven (2009) and for off-axis ob- 
servers by IVan Eerten et al.l d2010al) . However, due to the 
vast increase in numerical resolution provided by the boosted 
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FIG. 8. — Same as Fig. [7] now for v > v c . 

frame approach, this is the first time the post-break slopes 
have been determined from simulations where the jet break 
is fully resolved, and the current values can be considered 
quantitatively accurate. These slopes should be compared to 
observational data, such as the systematic study of Swift X- 
ray afterglows presented by iRacusin et al.l ([2009), that show 
a post-break slope for their sample of afterglows exhibiting 
'prominent jet breaks' that centers around a\ ~ 2. This dif- 
ference in slopes means that it is exceedingly difficult, at least 
for the Swift sample and at least for on-axis observers, to rec- 
oncile the data with a model of an initially top-hat blast wave 
decelerating into a constant medium. Even for off-axis ob- 
servers this is becoming problematic, although a number of 
caveats apply: the jet break might be simply post-poned be - 
yond what Swift can observe (Van Eerte n et al.l2010al l201 la), 
or only a fraction of an off-axis jet break is seen (see also Fig. 
[9j discussed below). Instead, the post-break slopes are more 
consistent with the values normally associated with high- 
lati tude emission ('reg i on I' of the 'canon ical light curve', 
see iZhang et al.ll2~006at IRacusin et all 120091) . without neces- 
sarily implying that these should be interpreted as such, since 
this interpretation would require extremely narrow jets (em- 
bedded in quasi-spherical outflow in order to get regions II- 
IV of the canonical light curve) and simulations of jets with 
Oo <C 0.05 rad., the smallest angle discussed in this paper. 
Possibly, Swift GRBs do not predominantly explode into a 




t„ bs (days) 



FIG. 9. — A comparison of sharp power law fits and synthetic light curves 
at observer angles o t, s = 0, 0.12, 0.2, rad. (top to bottom) for 9q = 0.2 rad. 
Plotted is the case v < v c . Other parameters are set as follows: p = 2.5, 
e e = 0.1, e B = 0.01, §jv = 1.0, z = 0, d 2 S = 1, n = 1 cm" 3 , E iso = 10 53 erg. For 
clarity of presentation, only half the data points of the synthetic light curves 
are plotted. Two of the three curves have been scaled by a factor ten, again 
for presentation purposes. 

homogeneous medium but in a different environment (e.g. 
stellar wind instead). Alternatively, the jet break might be 
hidden from view by an additional physical process, such 
as p rolonged injection of energy into the bla s t wave (see 
e.g. iNousek et al]|2006t Panaitesc u et al.l 120061 : [Zhang et al. 
2006b; Panaitescu & Vesu-anHojJ) - 

In Fig. [9] we show light curves and sharp power law fit 
results for Q \, s = 0.2 rad. These illustrate the extent to which 
sharp power law fits overlap with the data. In practice, it is 
not very difficult for off-axis observations to push the final 
turnover associated with the jet break out in time beyond the 
timespan typically covered by Swift (i.e. 10 days), especially 
once nonzero values for redshift z are considered, which can 
lead either to a missing jet break or a steepening that is far 
more shallow if detected at all because only the early part 
of the jet break is covered. An example of a jet break that 
is not fully detected is shown by the green dashed curve in 
Fig. [8] for 8 bs = 0q = 0.2 rad. However, in order to properly 
quantify these effects for e.g. Swift, an approach is required 
that includes not just synthetic light curves but also accurately 
models instrument biases a nd expected mea s urement errors , 
similar to the one taken by I Van Eerten et all d2010al 1201 lah . 
This will be the topic of a future study. 

After the onset of the jet break, the time evolution of the 
light curves in general does not follow a single power law evo- 
lution, as can be seen from Fig. [10] Given that the synthetic 
light curves consist of 85 (60) datapoints and were given arti- 
ficial measurement errors of ten percent, the reduced \ 2 val- 
ues for the various fits reported in tables [2] and [3] demonstrate 
that power law fit functions nevertheless fit the light curve sur- 
prisingly well. Even the v > v c fits for 9q = 0.05 rad. have a 
small reduced \ 2 value. The reason that these are neverthe- 
less noticeably higher than the other fits can be inferred from 
the late time behavior of the temporal indices for the narrow 
jet in Fig. [TOj Above the cooling break, the emission is dom- 
inated by a smaller region closer to the shock front than is the 
case below the cooling break. As a result, the observed flux 
above the cooling break at any given time consists of contri- 
butions from a smaller timespan in emission times. It will 
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FIG. 10. — Temporal index evolution for v < u c (top plot) and v > v c 
(bottom plot), for p = 2.5. 9 = 0.05, 0.1, 0.2 rad and 9 ohs = 0, O.60 o , &0< using 
the same colors and line styles as in Fig. [3] The top grey lines indicate the 
theoretical pre-break value. The narrow bottom grey lines indicate the ranges 
of sharp power law values found for the post-break slope, for all opening 
angles and observer angles except 9 B i, s > OAOq with 0q = 0.2 rad, where the 
temporal index did not reach a minimum during before 26 days. The thick 
bottom grey lines denote the middle of these ranges. 

therefore take less time for a change in the nature of the evo- 
lution of the blast wave to become noticeable than for flux be- 
low the cooling break, as illustrated by the comparison shown 
in Fig. [TT| What is seen for the v > v c curve at late times 
is the onset of the transition to the non-relativistic regime, a 
consequence of the fact that the smaller the opening angle, 
the smaller the total energy in the jets (with energy in jet and 
counterjet Ej « E, Sl 



5.2. Implications for the break times 

The evolution of the jet break time, as determined using 
a sharp power law fit, is shown in Fig. Q~2] If there were 
no lateral spreading at all, the jet break would be determined 
completely by the different edges becoming visible, and as 
a result the onset tjo and end tm of the jet break would be 
given by r fon cx (flr> - ^) &/3 a nd t h oc (0 o + obs )^ 3 respec- 
tively (Van Eerten et al. 2010a). For a jet observed on-edge, 
the nearest edge is visible already at r = 0, while the relative 
angle of the far edge is at its maximum distance of 29o. In 
reality, a steady spreading of the jet does occur. Also, the in- 
termediate light curve slope change at the onset of the break 



FIG. 1 1 . — A comparison between light curves for u > u c and v <v c and 
6q = 0.05 rad, o i, s = 0.0. Other parameters are set as follows: p = 2.5, e e = 
0.1, t B = 0.01, In = 1.0, Z = 0, d 2 S = 1, n = 1 cm 4 , E iso = 10 53 erg. For 
clarity of presentation, only half the data points of the synthetic light curves 
are plotted. 




FIG. 12. — Jet break times averaged over a range of p values for v < v c 
(top plot) and v > v c , as determined from sharp power law fits for different 
observer angles o t s and different jet opening angles 9q. The solid grey curves 
indicate r oc (^o + ^ofts) 8 ^ 3 - The dashed grey lines indicate break times for an 
on-axis observer, scaled from the on-axis break time to 05 for 9q = 0.05 rad., 
using r = 7-0.05 (0 O /O.O5) 8 / 3 . 
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is not as steep as the final slope change at the end of the break 
even for pure radial flow. These facts, together with the fact 
that the onset of the break is usually sufficiently early to be 
overwhelmed in light curve data (e.g. from Swift) by other 
early time features such as flares or plateaus, render it likely 
that in practice it is the end of the jet break rather than the on- 
set of the jet break that will be captured by a broken power law 
fit to the data. The relation between measured break time and 
jet opening angle will therefore lie closer to 77, oc (0o+$ofa) > 
than to Tb oc (#o) 8 > f° r general observer angle. Although the 
inferred jet breaks for the synthetic light curves do not fully 
reach this upper limit, Fig. Q~2] shows that, when the observer 
moves noticeably off-axis, they do trace this expected behav- 
ior at least for 6q = 0.05 rad. and 6*o = 0.1 rad. The jet break 
time as a function of observer angle is very noisy for the wide 
jet with 0q = 0.2 rad., mainly because in this case the jet break 
for observers far off-axis is not fully covered within the times- 
pan of 26 days. For small observer angles, when both onset 
and end of the break are still fairly close to each other, the two 
breaks have not yet fully separated and the turnover is still de- 
scribed by a single smooth break centered at r oc (0 o i, s ) 8 ' 3 , as 
indicated by Fig. [12] and the drop in a values for sB fits be- 
tween for increasing 8 t, s (as shown in tables[2]and[3]l. 

5.3. Implications for fit functions 

A comparison of the \ 2 fit results for the different fit func- 
tions shows that the performance of the different functions 
is comparable. Smooth power law fits of type sB by defini- 
tion outperform sharp power law fits, since the latter are a 
special case of the former, with a —> 00. Smooth power law 
fits of type sH perform poorly for on-axis observers, but of- 
ten slightly outperform the other fit functions for off-axis ob- 
servers, which is remarkable since fit function sB has more 
free parameters. Which fit function to use in practice will de- 
pend on the goal of the fit. If the goal is to obtain a fit as 
close to the data points and with as few parameters as pos- 
sible, sH is a good starting point. On the other hand, if the 
aim is to derive model parameters from the data for the type 
of model discussed in this paper, sB or even sharp power laws 
(PL) might be preferable, since especially for observers close 
to the axis, their ao values lie consistently closer to theoretical 
expectations (and model input for the synthetic light curves). 

6. SUMMARY AND DISCUSSION 

In this paper we present light curves for gamma-ray burst 
afterglows decelerating into a constant density circumburst 
medium. These light curves have been calculated from high- 
resolution AMR RHD simulations on a grid that is given a 
Lorentz boost in the direction of the jet, relative to the ori- 
gin of the explosion. The advantage of this approach is that 
the relative Lorentz factors in the outflow are reduced and 
Lorentz contraction of the shock front no longer presents a 
numerical resolution issue when blast wave deceleration at 
early times is calculated. The added complexity introduced 
by the loss of simultaneity across the moving grid relative to 
the rest frame of the burster can be dealt with by local inverse 
Lorentz transformations. A linear radiative transfer approach 
to synchrotron emission through the evolving fluid as repre- 
sented by a large number of data dumps from the sim ulation 
dVan Eerten & Wiiedl2009t IVan Eerten et alj|2010bl) is still 
possible, as has been presented in this paper. 

The dyna mics of narrow and ultra-relativi stic jets will be 
discussed in MacFad ven & van Eertenl (120121) . In the current 



study we focus on the radiation and the nature of the observed 
jet break. In a given asymptotic spectral regime, the shape of 
the light curve is completely determined by the scale-invariant 
evolution of the spectral breaks and the peak flux. The func- 
tions describing these evolutions are characteristic functions 
of observer angle o b S and initial jet half-opening angle #0 
only and can be scaled between different explosion energies 
and circumburst densities. Since they are also independent 
of synchrotron accelerated particle slope p, they can be used 
to generate light curves for arbitrary value of p. Generalized 
scaling relations for arbitrary circumburst density profiles (in- 
cluding ISM and stellar wind) are provided. 

The time evolutions of the spectral breaks and peak flux 
change directly following the jet break, although, thanks to 
the vast improvement in resolution, an earlier reported tem- 
porary post-break steepening of the cooling break v c is found 
to have been resolution-induced. Nevertheless, the tempo- 
ral behavior of v c for off-axis observers was found to be ex- 
tremely sensitive to small deviations from radial flow, even 
at early times and this is likely to leave an inprint in obser- 
vat ions, although any specific model (such as a structured 
iet.lMeszaros et al.l|1998l: iRossi et"aLll2002h iKumar & Granotl 
12003b iGranotl 120051) prediction might be hard to disentangle 
from the effe cts of changes in the synchrotron emission pro- 
cess (see e.g. Filgas et aC 20Tlb. 

The shape of the jet break is systematically surveyed for 
jet opening angles 6q = 0.05, 0.1, 0.2 rad., observer angles 
ranging from on the jet axis to on-edge and p values rang- 
ing from 2.01 to 3.0. Pre-break temporal indices are found to 
be in good agreement with theoretical expectations for purely 
radial flow. This is partially a consistency check on the com- 
puter code, since purely radial BM flow was used to set up 
the initial conditions of the simulations. On the other hand, 
the simulations were started from an ultra-relavistic on-axis 
Lorentz factor 70 =100 and minor deviations from radial flow 
will therefore have occurred well before the jet break. 

For the cases considered, post-break temporal indices are 
generally far steeper than theoretically expected for a quickly 
expanding jet. This doe s not imply exponential jet ex- 
pansion actually occurred (I Van Eerten & MacFad yenll20 1 2bl 
MacFadve n & van Eertenl 120121) as demonstrated by the de- 
pendency of the jet break shape on the observer angle, but 
represents the combined effect of steady expansion and the 
edges of the outflow becoming visible. The difference in 
slopes between the s ynthetic light curve s and those reported 
for the Swift sample ( Ra cusin et al.ll200"9l) means that it is ex- 
ceedingly difficult, at least for the Swift sample and at least 
for on-axis observers, to reconcile the data with the model 
of an initially top-hat blast wave decelerating into a constant 
density medium. Even for off-axis observers this is becom- 
ing problematic, although a number of caveats apply: the jet 
break might be simply post-poned b eyond what Swift can ob- 
serve (I Van Eerten et a l. 2010a, 201 la), or only a fraction of an 
off-axis jet break is seen. Sharp power law fits confirm that the 
jet break time is sensitive to the observer angle and increases 
significantly as the observer moves off-axis, which has im- 
plications for the interpretation of afterglow data and inferred 
energy of the explosion (which will be overestima ted when an 
on-axi s observer is assumed, as discussed in I Van Eerten et all 
I2010al) . 

This discrepancy between light curve slopes from ISM sim- 
ulations and Swift (or other instrument) data can in theory be 
explained by assuming that afterglow blast waves decelerate 
instead into a stellar wind environment shaped by the progen- 
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itor star. The most likely scenario then is one where a jet in 
a stellar wind environment is viewed almost on-edge, given a 
random orientation of the jet. The jet-b reak is generally less 
steep for a stellar wind environment (|Kumar & Pana itescu 
l2000tlGranotll2007HDe Colle et al.ll2012bl) . A further compli- 
cation is added by the fact that GRB progenitor stars are not 
expected to exist in complete isolation, and the stellar wind 
environment of the sta r is likely to be shaped by multiple 
colliding stellar winds (Mimica & Gia nniosll201 ll) . Full re- 
sults for the stellar wind case computed from a boosted frame 
will be presented in a follow-up study. Alternatively, the 
explosion does occur in a homogeneous medium but the jet 
break is hidden from view by an additional physical process, 
such a s prolonged injection of energy into the blast wave (see 
e.g. iNousek et alj|2006t i Panaitesc u et al J 120061: IZhang et all 
2006bl iPanaitescu & Vestrandll2012l) . ~ 

Different power law fit functions have been used in the 
literature to describe jet breaks. Compari ng sharp power 
laws, smoothly connected power laws (\sB\ Beuerman n et alj 
fl999l) and power law transitio ns including an exponential 
term ('si/'. lHarrison et al. 1999), we find that all descriptions 
provide good fits to synthetic light curves, although type sH 
underperforms for on-axis observers and often outperforms 
the other types for off-axis observers. Nevertheless, type sB 
fit functions and sharp power laws yield pre-break results that 
are the easiest to interpret in terms of the underlying model. 

The simulation data, light curves and characteristic func- 
tions (i.e. the scale-invariant time behavior of peak flux and 
spectral breaks) in this work will be used to improve the ac- 
cur acy of simulation based data fitting methods such as BOX - 
FIT jVan Eerten et alj2012l) . From the characteristic functions 
light curves for each asymptotic spectral regime can be repro- 
duced directly, as has been done in this paper and f ollowing 
the approach from lVan Eerten & MacFa dyenl (120 1 2al) . For the 



full spectrum, a heuristic description of t he sharpness of spec- 
tral transitions is re quired as well (see Gra not & Saril 120021 : 
iLeventis et al.l 120121 for an example of this approach in the 
spherical case). Alternatively, the simulation output for the 
blast wave dynamics can be processed using the methods em- 
ployed for BOXFIT, albeit with the extra step of transforming 
to the lab frame. This has the advantage that radiative trans- 
fer equations can subsequently be performed very quickly 
and that no heuristic description of the spectral transitions is 
needed. 

As stated earlier, the steepness of the post-break slopes 
poses a challenge for the Swift sample. A true test of the sever- 
ity of this issue is to compare observational data and synthetic 
light curves systematically using one of the simulation-based 
fit approaches described above. This will be the topic of future 
work. We note that the one afterglow that has already been 
fitted using the BOXFIT approach, GRB 990510, has a steep 
post-break temporal slope com pared to t hose in the Swift sam- 
ple (F oc f" 2 40 , according to Stanek et afl ll999h . which helps 
to explain how it was possible to obtain a good fit using the 
ISM model for that particular burst. 

All light curves and spectral break and 
peak flux evolution functions from this work 
will be made publicly available on-line at 
|http : //cosmo . nyu . edu/af terglowlibraryj 
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APPENDIX 

LIGHT CURVES FROM A BOOSTED FRAME 

As in previous work (I Van Eerten et alJfeOlOb a), the radiative transfer equation is solved in the burster frame (the "lab" frame) 
simultaneously for a large number of rays through the evolving fluid. Because the simulation grid is itself moving with a fixed 
Lorentz factor, the simulation frame is no longer equal to the lab frame. The consequence of this is that additional Lorentz 
transformations will be necessary going from simulation to burster frame, not only to boost the fluid quantities, but also to take 
into account the loss of burster frame simultaneity across a single snapshot. As a consequence of the latter, the contributions from 
a single snapshot to the emission and absorption coefficients of the rays, for a given observer time and angle, no longer lie on a 
flat intersecting plane (previously labeled 'equidistant surface' or 'EDS'), but on a curved surface. Denoting observer time t„b s , 
burster time t and simulation grid time f ', we have for each ray on each snapshot the constraint 

t obs = t(t')-R(t'), (Al) 

where R the distance traveled by the ray in the burster frame parallel to the line of sight. R is equal to zero when the ray crosses 
the EDS plane centered on the burster frame origin and oriented perpendicular to the line of sight (i.e. defined such that light 
emitted from the origin at t = will arive at observer time t t, s = 0). This plane, which we label 'EDS0', is defined in the burster 
frame and therefore still flat.). For any given ray, the relevant coordinate for a given snapshot is 

q = q E + Ru E = q E + (t-t oh ,)u El (A2) 

where q E the coordinates of the point where the ray crosses EDS0 and u E a unit vector pointing along the ray to the observer. 
Writing the vector components of the previous equation explicitly, we get 

/ 9x \ ( lEx \ ( sin9 ohs \ 

% = lEy +{t-t„ bs ) . (A3) 

\ 1z J \ 9Ez J V COS& obs J 

Before the radiative transfer calculations are performed, we pre-process the snapshot files to store the local fluid states in terms 
of burster frame coordinates (q,t). The relevant Lorentz boost equations for a boost of factor 7s and velocity j3$ along z, the 
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direction of the jet, are 



q' z = j(q z -l3st). 



(A4) 



Combining these with equation IA3I allows us to determine which fluid cell to probe for a given ray (determined by its qn 
coordinates) and given snapshot (determined by its simulation frame time t '), leading to: 



Ct' COS 6 bs 



ct obs cos 6 obs 



1 - fis cos 6 obs 7 S ( 1 - fa cos 6 obs ) 1 - /3s cos obs '' 
Is c 

q y = qEy 

q x = qEx + c(t - t obs ) sin 9 obs . 
The distance dR traveled by each ray between two snapshots that are df ' apart is given by 

cdt' 



dR = cdt = 



7s(l-^ s cos6»„ fcs )' 



(A5) 



(A6) 



The local emission and absorption coefficients are a function of comoving fluid number density n, comoving fluid energy density 
e and burster frame fluid velocity v. The comoving quantities are provided directly by the fluid simulation, since they are 
independent of the grid velocity. The velocity and Lorentz factor in the burster frame are calculated during the pre-processing of 
the grid snapshots according to the standard relativistic velocity addition rules. 
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